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ABSTRACT 


This thesis investigates the pressure drag coefficient in 
the transonic regime over an axi-symmetric body, witha set of 
unique contour surfaces developed in a previous thesis. The 
contour surfaces were obtained by an exact solution of the 
small perturbation transonic equation, using the guidelines 
and tools developed at NPS. In this work, Computational Fluid 
Dynamics (CFD) was not only used to compute the afterbody 
contour surface, but also to investigate a conical afterbody 
and complete bodies, which are composed of an arbitrary 
forebody (ellipsoid) and variable afterbody (contour and 
conical). Euler as well as Navier-Stokes flow-solvers were 
applied to the geometries of interest, giving Mach-number 
contours for viscous and inviscid flow, pressure drag 
coefficient magnitude, and depicting shock wave locations. On 
the basis of these results, it can be verified that our 
contour surface afterbodies will decrease by 15% the peak of 
the pressure drag coefficient (Cy) versus Mach number curves 
in the transonic regime. These results can be used to design 
low pressure drag surfaces for such shapes as missiles, 


projectiles and aircraft engine nacelles. 
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L. INTRODUCTION 


The flight of an object over a wide range of speeds has a 
Critical transition zone where both subsonic and supersonic 
myoes of flow exist. This speed regime is referred to as the 
transonic range. The critical aerodynamic behavior occurs in 
the range 0.8 < M < 1.2 depending on the object, where the 
aerodynamic coefficients have been found to change by as much 
as 100%. The behavior of the aerodynamic force components is 
usually characterized by a rapid increase in the coefficients 
followed by a sharp drop, in other words, a peak value for the 
coefficients arises in this regime [(Ref.1 and 2]. 

One of the aerodynamic components, aerodynamic drag, 
represents a significant adverse force on all flying objects 
such as aircraft, missiles and projectiles. A high drag force 
reduces the craft’s range capability or equivalently requires 
more energy to achieve a certain range. Any effort to reduce 
the drag coefficient in the design process must concentrate on 
reducing the wake and pressure drag (inclusive of wave drag) 
@emernibutions to the total drag. 

A traditional approach to investigate aerodynamic 
characteristics is based on wind tunnel data and actual flight 
testing as a source for improvement of configurations to get 
the best results. Unfortunately, both wind tunnel and flight 


tests are considerably expensive and time consuming, 


particularly in the transonic flow regime. Because of the 
great lateral influence of profiles at transonic speeds (shock 
reflection of the wall bounces back to the model), models must 
be made extremely small compared with wind tunnel dimensions 
and this introduces great experimental difficulties. 
Consequently, the data obtained at transonic speeds is 
considerably less reliable than at either subsonic or 
Supersonic speeds. In addition, from a mathematical point of 
view, even the two dimensional small perturbation potential 
equation for transonic flow retains one non-linear term which 
is essential for non-divergent solutions at Mach one, but this 
non-linear partial differential equation has proven to be 
Gifficult to solve. Such inherent difficulties, coupled with 
the presence of shocks in the flow which cause boundary layer 
separation, have resulted in the creation of many approximate 
methods of solution which are employed in the design of 
transonic airfoils and the like [{Ref.3, 4 and 5]. On the other 
hand, the use of numerical simulation known as Computational 
Fluid Dynamics (CFD) to predict aerodynamic characteristics 
greatly increases possibilities to improve design optimization 
at relatively low cost and allows for ease of design changes. 

Finally, using the latest capabilities of Euler as well as 
Navier-Stokes flow-solvers, it has been possible to compute 
the flow over axi-symmetric bodies with various contours in 


the transonic regime. 


II TRANSONIC FLOW 


Transonic flows are characterized by the simultaneous 
presence within the flow field of both subsonic and supersonic 
regions. The properties of transonic flows can be derived from 
the general equations of gas dynamics, namely, the equations 
of state, continuity, momentum and energy. The following 
derivation based on assumptions that the flow is steady, 
irrotational and isentropic with no energy transfer, no body 


fereces and no shear stresses (inviscid flow). 


1. Small Perturbation Theory 
Starting with the equation of motion for. steady, 


isentropic, inviscid flow in the index notation form [Ref.4]: 


(1) 





Small-Perturbation Theory gives the velocity field as 


a Ue, oe U, = W (2) 


where U denotes the free stream velocity in the x direction 
and u, v, w are called perturbation velocity components in the 


direction x, y, z respectively. 


Writing eq.(1) out in full, and substituting the Wetoeme 
field eq.(2), gives the equation in terms of perturbation 


velocities. 


Gu, OV , OW) . gE i a 
= ee Oy i 5a Tau) Ox 4 Oy ‘ OZ 
(3) 
Ou Ov OV Ow Ow Ou 
+ (Utu) oe + =) & re ~ = + w(Ueu) (Se - =a 


From the energy equation for a perfect gas, the speed of sound 


(a) can be expressed in terms of the perturbation velocities. 


Cea A ce, Bs (4) 
2 iy a5 2 (yee!) 
On 
a? =a - VY) (2uu + u’ + ve + we) (5) 
Substituting eq.(5) into eq.(3), dividing by a ema 


rearranging the terms, gives the full exact equation in terms 
of perturbation velocities and free stream Mach number. This 
equation contains linear terms on the left-hand side, but on 


the right-hand side the terms are nonlinear. 
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Bu. 2u) . My. uy du, du) , vl du, 2v) 
OZ "Ox U2\ oy dz 








If the perturbation velocities are small (u/U,v/U,w/U << 1), 
eq.(6) can be simplified by neglecting the terms containing 
squares of the perturbation velocities on the right hand side, 


yielding 


2, OU OV Ow 
SS ae 
2 u Ou oe ly OVE. 
Mey) oe + Mo (y 1) =— 3) (7) 


ee 


2 V/ Ou Chia 2 W/ OU Ow 
| OZ xx | 


For further simplification, in eq.(7) all the terms on the 
right hand side can be neglected, in comparison to those on 
the left hand side. This gives the linear equation, which 
contains only perturbation velocities and is valid only for 


subsonic and supersonic flow. 


Clee a Oa (8) 


For transonic flow, where M, > 1, the coefficient of du/dx on 
the left hand side becomes very small, but it is not correct 
to neglect the first term on the right hand side of eq. (7). 

Therefore, the governing equation of transonic flow in term of 


perturbation velocities is as follows: 


20 OV , OW _ 42 


> . a Ree ule: 
(1-Me) Set et ae = Ma (yl) 5 5e (9) 


For irrotational flow, a perturbation velocity potential @¢ 


exists, 
dy dy dy 
- ~ = 10 
US Fx er ener oe 
Substitution of eq.(10) into eq.(9), gives the governing 


equation for transonic flow in terms of the velocity 


potential. 


2) A“ " 0“ "7 0° a Me (+1) dy d*0 ‘Soe 
ax? dy? dz? U OX 9x? 


For bodies of revolution, it is convenient to use cylindrical 
coordinates (x,r,@) where x is aligned with the body axis. The 


Metlocity Components. corresponding to (x;r,@?) are u,, u, and 


1% 
Maa respectively. The velocity potential in the cylindrical 


coordinates, 


(1:2) 


Transforming Cartesian coordinate eq.(9) into cylindrical 


coordinate, gives the governing equation of transonic flow. 


ee ee cE 


Lod, 1 0G _ Bett 96 9° (13) 
Q@x2 dr? ror y2 9g? 0X 952 








For axially symmetric flow, where the conditions are the same 
in every meridian plane, there is no variation with @, so the 
small perturbation, non-linear, axi-symmetric transonic 


potential equation can be written as follows 


A°h _ Me(y+1) g a6 d 


i sxeguy “¢ (14) 
Q@r2 ror °" ax? U OX 9x2 





Or rewriting in shorthand notation, 


72 
Meas 1) 
bag bn + (1M Ge (15) 


Reference 6 introduces a modified potential equation for axi- 


symmetric flows 


1 0 2 
Saye (br) + (1-MS) bax = bb xx (16) 


where the modified velocity potentials are 


be = ME (yt), bp = Ma (y*1) (17) 


2. An exact solution for axi-symmetric, transonic flow. 
Solutions to the modified transonic eq.(16) have been 
given by Biblarz [Ref.6 and 7] by using the separation of 


variables approach with a potential function ¢(x,r) of the 


EO 


6(x,r) = &(x)n(r) + (1-MS)x (18) 


Substituting the above ¢@ function in the modified transonic 
eq.(16) results in two ordinary, second order, non-linear 


differential equations 





dé d*é | _ 
GEE NE 0 (19) 
and 
d*n eal An? = 0 (20) 


where A is the separation constant. 
The solution to the first differential eq. (19) is obtained 


by multiplying both sides dé/dx, 


dé |dt d*t| _ dé z (21) 
dx |de qa |~ ae\“*) = ° 
or 
ah ey _ Az2 |= 9 (22) 
dx} 3\dax Z 
Thus 
dé 3 ; 
% ee (23) 
Se ( gM? +a] 
where @w@ iS a constant given by Ref.6é. 
Rearranging eq. (23) 
dx = —_ 
‘e (24) 
[ SA2? pe al 3 


and integrating eq. (24) gives, 


dé 


pe 


2 Gee 
2 (25) 


where @ and X, are integration constants. 


The solution to the second non-linear ordinary 


differential equation (20), is obtained by an outer expansion 


method. [Ref.8] 





n(x) = + (1-MZ) £,(r) + (1-MZ)?2£5(r) +... (26) 


Where (1-M,°) represents a small parameter and the first term 
is the purely sonic solution. 


By taking first and second derivatives eq.(26) and 


subtituting them into eq. (20) yields, 


(1-M3)? 


n(x) ==, +|1-M2 | adr ‘VB +2) 4 22,272 (vo+2) , (C27) 


28+8/8 


Al-Hashel [Ref.9] reported on eq.(26) and eq. (27), 
implemented the boundary conditions with the consStaneeua ae 


a and \, gave the final results in the new variables, 


(28) 





plVB+2) = -(VF+2) yn 2 (29) 


_ i 1 
zac? 43h)? (30) 
* 2 
a n A 
a (a ,) 0-4142 (31) 
then 
a dé 
J i (32) 
where X)=0, and 
1 2 OT +2 (1-Mz)? O+ (33) 
ee | a Meae yee), Ae ee) 
x2 50.63 


Equation (32) has been numerically integrated and plotted on 
figure (1) as (xX) versus X, and eq. (33) will be evaluated and 
mlOeted in figure (2) as #(r) versus fF for M, = 1.05, 1.1 and 
1.2. A "patching" technique discussed in Ref.9 has been used 


here. 
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III. PRESSURE COEFFICIENT AND BOUNDARY SURFACES 


1. Pressure Coefficient 


Liepmann and Roshko [Ref.4] define the pressure coefficient 


as; 


‘a Vee) & 2 1D . 
> P : (34) 


From the isentropic relation, we have pressure ratio in terms 


of Mach number and after substituting into eq. (34) yields 
5 W/(y-2) 
Pe 2 | 2s (35) 
pO 3 =: 


yMe ||2+(y-1) M7 


Introducing M,”? = U*/a,*, M* = u?*/a* and using energy equation 


(recall eq.4 and eq.5), the pressure coefficient can be 


expressed in terms of the perturbation velocities, 














FED, 
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(36) 
J itys2) 

_ 2 _¥71y? 2u, ut+ve+w? a 

y) 2 = ae 72 
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eZ 


Using the binomial expansion on the expression inside the 
square brackets in eq. (36), we obtain the pressure coefficient 


in the form 


2uU 


2 Pa 2 
gine Se ee meas (37) 





Por axi-symmetric £Elow, 1n cylindrical coordinates where u=u, 





and (v*+w“) = Clee , substituting into eq. (37) yields 
Ps 2 
Ce Bes: 7 (awl 22) ated (38) 








The linearized pressure coefficient approximation for axi- 


symmetric flow turns out to be 


C, = -— (39) 


Recall the modified axial velocity potential, eq. (17), 


u 
thus 
u p 
i 7 aa (41) 
ao fags) 


Substitute eq. (41) into eq.(39), yields 


-2¢ 
Cp = ———— (42) 
Nis Wy 2ab) 


The derivative of the potential function eq. (18) with respect 


to x becomes 
(43) 


Ref.6 introduces relation of the constants C,, a, a and Xd in 


expression, 


a = +C,|1-me|1-7574 (44) 


and 


(45) 
es LOE 


Ci -2 
ee eel) 
r \0.7574 





Rewriting eq.(23) and substitute the constant a eq. (44) 


becomes 


1 
dé 2 iw 3h 2 3 4 
= : Shy? + Cy] (1-2) |2-7574]3 (46) 
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Recall the new variable eq. (28) 


reer 
rl 


_ (3 


(47) 
es 








Substituting eq.(47) into eq.(46) and factoring C, out, yields 


1 
a C, [E2 + | (1-M?) p2-7574/3 (48) 
Recall and arrange new variable eq. (31) 


0.4142 
Seal mee (49) 


Rewriting eq.(43) the modified potential function 


inserting eq. (48) and eq. (49), 


and 


becomes 


We) 


1 
0.4142 : bs 
ine (ah) eh E2+| (a-m2) |2-7574]3 - (1-m2) (50) 


or 


1. 
qi-4142¢ * 


1 
ae 2 z (S51) 
\ 0.5858 d | g*+| (1-Mg) pom + (1M) 


o 


Px = 


iS 


Substituting the expression of constants eq. (45) into eq. (51) 


1 
3 2 (52) 


j- s= 02208 i £2 + | (1-M2) hove + (1-M3) 


Finally, rewriting the pressure coefficient C, as 


a 
C. = _-2_b aoneaete (1-M2) 7874] F a’ C53) 
Mo tye) 


2. Boundary Surfaces 

For an inviscid flow, the condition to be applied at the 
surface of a solid boundary is that the direction of the flow 
velocity vector must be tangent to the solid surface. In other 
words, the velocity vector is everywhere at right angles to 
the normal to the solid boundary [Ref.10]. In addition the 
boundary condition requires that the gradient of potential ¢ 
vanish far ahead of the body. In terms of perturbation 


velocities this boundary condition becomes 


u 
ee 4 
7a (54) 


e 
AX 


a 
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Recall the modified velocity perturbation potential eq. (17) 


u 
br = MS (y+1) (55) 
Ue 
thus 
u p 
= == (56) 
co MS aye) 
substituting eq. (56) into eq. (54), 
obe p 
(S) = ee (57) 
Mee yd) 


Taking the differential of eq.(18) with respect to r gives 


d 
o-=¢& = (58) 


Expressing eq. (58) in terms of new variables, recall eq. (28), 


eq. (29) and eq. (31) 


ao an 
o,- = 0.0848 ~ =~ (59) 
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Recalling eq. (33) and taking the derivative with respect to f 


then substituting eq.(59) into eq.(57), yields 


ee = DSUSEeEe Bae Boy ja-Mo| £1-828440.1512 (1-M5) ee 
aX} Mo(y+1) | F? 


(60) 


Replacing the left hand side of eq.(60) with the new variables 


eq. (29) and eq(30), then rearranging it becomes 
OF = 0.0526 F | “8.2 g2g4 |1-M2| £1-8284+0,1512 (1-M2)? £6-657 
aK Me (y+1) | £? 


(61) 


For further arrangement, by separating variables eq. (61), 


yields 


ar 
8 


= 3 


2.8284{1-M2| £1-8284-9 1512 (1-M2)2 F&-657 
(62) 
-0 .0326 


g dk 
Ms (y+1) 


Then, integrating both sides of eq.(62) becomes 


Ite! 


; drt 
| [Ea oases m] 


5 -2,8284|1-m2|F1-8284 -9,1512 (1-M2)? F6-657 
aL 


(63) 
0.0528 [ & az 
MOV tL) 
where 
ee OE 
0 Ja -m2 |°-2072 (64) 


Al-hashel [Ref.10] has developed and computed eq. (63) 
using numerical integration, to determine the boundary 
surfaces in dimensional and non-dimensional (normalized) form 
tome M = 1205, 1.10. and 1.20 as depicted in Figure (3) and 
Figure (4). Based on these calculations, the geometric grid of 
the afterbody for Mach number 1.10 and 1.20 are developed for 
further study with Computational Fluid Dynamic (CFD). 

This thesis research also examines conical afterbodies as 
a solid afterbody boundary surface with base diameter ratio 
(rye) OL O50 and 0-75 and conical turning angle (8) of 26.6 
and 14.0 degree respectively [Ref.11 and 12]. Then, for 
further investigation, the complete bodies as a solid boundary 
surface are generated, with a kind of forebody (ellipsoid) 
joined with the contour surface afterbodies as well as the 


conical afterbodies. 
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IV. DRAG 


Drag is one of the aerodynamic force components parallel 
to free stream velocity (U,). It represents a significant 
adverse force on all flying objects. Basically, the drag force 
is divided into two categories, the drag caused by forces 
acting normal to the boundary surface which is called pressure 
drag (inclusive forebody, base and wave drag) and that arising 
from the tangential forces acting on the surface, by virtue of 
viscosity which is called viscous drag or skin friction. 
Figure (5) shows the components of the drag as function of 
Mach number and the methods used to compute body drag in four 


Mach number regimes [Ref.5]. 
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Figure 5. Components of the Drag [Ref.5] 
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ie ‘okin Friction 

Skin friction is the result of shearing force within the 
boundary layer of a viscous flow, acting tangentially to a 
surface in motion relative to the fluid. The amount of viscous 
resistance depends upon whether the flow is laminar or 
turbulent [Ref.5]. Krasnov [Ref.13] introduces the laminar and 
turbulent skin friction for flat plates, based on the boundary 
layer theory, for compressible flow. These formulas also valid 


for bodies of revolution with infinite length. 


Laminar flow, Re < 10° 


ae 
CA = ——— (1 + 0.03MS) 3 (65) 


(Re) 2 


Turbulent flow, Re > ape 


ae See 2 Meg raereO doa) ep = (66) 


(1+0.1M2) 


where Re is the Reynold number 


fais 


2. Pressure Drag 
Liepmann and Roshko [Ref.4] introduce the pressure drag 
formula (inclusive wave drag and base drag) for axi-symmetric 


bodies as, 


L 
D={ pads pO kee) 
Q 





(67) 
iE 
=f (p-p,) dS + (p,-p,) S(L) 
0 
or, in dimensionless form, 
1B, 
Ca 
my mS IE) 
(68) 
1 |. ds 
= (Go Sees (ob oe 
S(L) J Pqxo* * “pp 
where S(x) = mr* is the cross-sectional area of the body at x. 
= = dr 
GS'(30) = 2 Tiseia = ite (69) 
ax 
or 
— = amr Se (70) 


ns 





The first term in square bracket eq.(68) can be solved by 


inserting eq. (53), eq.(60) and eq.(70) yielding, 


a 
a eer [E2+| (1-m2) | 12-7574) Fea’ | 
Moy) 


——— Ab 
=P 
sary 
iT 
—— bh 


(71) 


(277) UGE = +2.8284|1-M2| #2-828449 1512 (1-2) 2#6- 657 hae 
Me (y*l))) |= 


Recalling new variables eq.(28), eq. (29), 


inserting B* = 1-M,? into eq.(71) yields, 


ds a 
ea 
1 
x FI z 
5 ie 52 1302C 
2 |p a20a9tF+1 93-5248) ) F482 nee 
(1-87) (y+1) V0. 7072 g0-2071 


(72) 


eee 2 8284 | B2| #1-828440 151284 #6 -657 lax 
(io © ea nag gnc a a 


where C,, a and A are constants. 
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The second term in square bracket eq.(68) is the contribution 
to the drag by the base pressure p,. Ref.4 claims, tEhat@eene 
values of the base pressure coefficient C,, must be obtained 
experimentally. However Krasnov [Ref.13] introduces’ the 


boattail drag coefficient for a conical afterbody as, 


Cpp = 0.002 [0.8 + — |61-7(1-Sp) (73) 
M. 
where 
6 = turning angle Of sence wcontcal bed, 


Sp 


ratio of the base area to the mid-section area 


In this work, the skin friction coefficient as well as the 
pressure drag coefficient of the geometries of interest, at 
specified Mach number, can be obtained directly from the 
results of CFD calculations. The Euler solution will give the 


pressure drag and the Navier-Stokes solution will give both. 
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V. COMPUTATIONAL FLUID DYNAMICS 


The rapid advancement in the speed of computers and their 
enormous memory size has led to the emergence of the field of 
computational fluid dynamics (CFD). This branch of fluid 
dynamics complements itsS experimental and theoretical 
branches, by providing an alternate, cost-effective means of 
Simulating real flow. It also offers the only means of 
examining theoretical advances for conditions unavailable 
experimentally. As a model based method, CFD can provide the 
convenience of being able to switch off specific terms in the 
governing equations [Ref.14], so as to assist the researcher 
in understanding the contributions of various physical 
Eeker@uae 

In this work, CFD was used to compute the axi-symmetric 
flow over the afterbody geometry of models only (the boundary 
surfaces obtained by the small perturbation method [Ref.9] and 
conical afterbody) and over complete body models which are 


composed of forebody (ellipsoid) and afterbody. 


1. Grid Generation 

The computer programs GRAPE ([Ref.15] and GRIDGEN2D 
[Ref.16] are tools used to generate two-dimensional structured 
grids about airfoils and other shapes by the use of algebraic 


or Poisson differential equation solvers. GRAPE was used for 
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the geometry with smooth contour surfaces, while GRIDGEN2D was 
used for the geometry with a conical afterbody with non-smooth 
points. Outer and inner boundaries were specified as the C- 
type grid for afterbody models only, while the O-type grid was 
for the complete body models, where both type of grids treat 
the surface of the body as the inner boundary. The important 
characteristics in a grid generation technique are the ability 
to specify the spacing between mesh points at the boundary, in 
the direction normal to the boundary, and the control of the 
angles with which mesh lines intersect the boundaries which is 
known as orthogonality [Ref.15]. 

Figure (6) is a typical output of program GRAPE for an 
afterbody with small perturbation solution contour for Mach 
1.10 (SPS_1.1), C-type grid with the grid size of 1i5x3e0 ene 
figure (7) is the complete body with the grid size of 152x60 
(O-type grid). 

To develop a three-dimensional grid from two-dimensional 
grid (output program GRAPE), the FORTRAN code called HALF.F 
(Appendix) was used to write out the half of the 2D grid and 
convert them into 3D plane-grid. Then, the FORTRAN code 
ROTATEGR.F (Appendix) was used to rotate the 3D plane-grid for 
11 planes plus 2 more to generate a 3D volume grid of half of 
an axi-symmetric body surface as shown in figure (8) and 
figure (9) with” grid” size “of 58x13x60 and 77xeeaa0 
respectively. These grids are ready for further processing 


with a flow-solver. 


ZA 


Program GRIDGEN2D gives the pattern of how to build 2D- 
grid. It consists of four subfaces which are treated as 
boundaries. Grid spacing is determined by setting up the 
Gistribution points of each pair subfaces. It uses equal 
Spacing along the body surface, while the direction normal to 
the body surface uses geometric spacing with specified width 
in the beginning. Figure (10) and Figure (11) are typical 
outputs of GRIDGEN2D for an afterbody and for a complete body 
with the grid size of 58x60 and 77x60 respectively. FORTRAN 
code D2D3.F (Appendix) is used to convert 2D grid output from 
GRIDGEN2D into 3D plane-grid. Then we apply the FORTRAN code 
ROTATEGR.F to generate a 3D volume grid of half of an axi- 
symmetric body as depicted in Figure (12) and Figure (13) with 
miemeerd Size of 5S8x13x60 and 77x13x60. 

In this research, we also attempted to develop a fine grid 
for complete bodies (O-type grid), where the radius of outer 
boundary is set up to be five times of the body’s length. The 
steps to obtain 3D volume grid are similar to those described 
above with the final grid size of 77x9x120 and the final 
geometry of a quarter of an axi-symmetric body instead of a 
half body. Figures (14) and (15) are typical fine grids of the 
complete body with a conical afterbody and a complete body 


with contour surface afterbody. 
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2. Flow-solver 

The OVERFLOW program [Ref.17] was developed by NASA -Ames 
Research Center. It uses either 3-D Euler or Navier-Stokes 
flow-solvers for inviscid/viscous flow, by setting True or 
False the parameter VISINP (viscosity input) in the input file 
(overflow.in). Before flow-solver code (OVERFLOW) is applied 
on the grid file, a formatted 3D grid file named "grid.for" 
must be converted into an unformatted grid file named 
"grid.in", by using FORTRAN code called READX.F (Appendix). 
Then, the NAMELIST input file parameter specification must be 
written for running OVERFLOW, it is called "overflow.in". The 
input parameter consists of the number of iterations, 
timesteps, calculation methods, smoothing, type of flow and 
boundary conditions for each grid. The value of angle of 
attack (ALPHA) depends on the orientation of the grid in the 
coordinate system. In this case, ALPHA is 180° (flow comes 
from the x-positive to the x-negative direction). The boundary 
conditions depend on the geometry corresponding to the final 
3D volume grid. Both input parameters, the angle of attack and 
the boundary conditions for each geometric shape, are 
tabulated in Table 1. 

The file overflow.in (Appendix) is a typical input 
parameter specification of the axi-symmetric body with a 
conical afterbody using the grid size of 77x13x60. This file 
input uses NSTEP=100, Mach number=0.95, ALPHA=180°. The 


calculation method depends on a central difference Euler term 
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in J, K, L and ARC3D diagonal factorization and dissipation 
scheme (IRHS=0, ILHS=2, IDISS=2). For the first attempt a time 
step DT=0.1, ITIME=1 and CFLMIN=5.0 was utilized. The boundary 
condition consists of 6 boundaries (NBC=6), with the type of 
BC IBTYP=30 (inflow) in the J positive (1), IBTYP=15 (Axis K 
round) in J negative (-1) direction, IBTYP=12 (symmetry in Y) 
in the K poSitive (2) and K negative (-2) direction, IBTYP=1 
(Inviscid adiabatic wall) in the L positive (3) direction and 
IBTYP=32 (Supersonic /subsonic inflow/outflow) in the L 
negative (-3) direction. 

The OVERFLOW program gives output files such as ovr.out, 
G@meewe, resid.out and Fomo.out. To verify that a calculation 
is appropriate or converged, it can be traced by looking at 
the plot of residual history (resid.out). A convergence 
criterion was defined as the reduction in residuals by two 
orders of magnitude. If the first run (NSTEP=100) does not 
fulfill the convergence criterion, one may perform a restart 
by further running OVERFLOW and changing the parameter 
RESTART=.T. in the input file (overflow.in) and by copying 
q.save into q.restart. Repeating the above steps until all 
convergence criteria are met. If in the output file q.bomb 
appears, one tries another run by changing parameter time step 
DT and CFLMIN until the criteria are met. Then, one converts 
the output file q.save (unformatted) into formatted file named 
q.form by using FORTRAN code READ1.F (Appendix). 


Finally, the last step is to run PLOT3D program with the 


Zo 


source formatted files x=grid:-for and q=q.form. “Eyes 
functions on the PLOT3D program, many plots such as pressure, 


velocity, Mach number, vorticity, etc., are obtained. 


3. Results 

The results of CFD programs are grouped into the 
corresponding geometric shape, namely, afterbody only and 
complete body. The Euler flow-solver (inviscid flow) was 
applied to all axi-symmetric bodies, except for the afterbody 
models only, where both Euler and Navier-Stokes (viscous flow) 
were applied. Most of the calculations converged in 500 
iterations, meaning that the residual history achieved a two 
order of magnitude drop. 

Results can be analyzed by plotting the q.form file, 
output file from the OVERFLOW program, using the proper 
function in PLOT3D program (Mach number). By interpreting the 
Mach number contour surrounding the body surfaces, one can be 
determine the characteristics of the flow field. In addition, 
the drag coefficient C,, can be obtained in the file femerene 
output from OVERFLOW. Hence, in the sequences of Mach number, 
one can describe the significant flow characteristic of each 


geometric shape. 


a. Afterbodies 
An Euler as well as a Navier-Stokes flow-solver was 


applied to these afterbody models. The approaching free stream 
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Mach number (ranging from 1.05 to 1.50) starts from the mid- 
section of the afterbodies, then the flow follows along the 
afterbody surface until it reaches a maximum local Mach number 
as tabulated in Table 2. The maximum local Mach number for 
viscous flow is always lower compared to in the inviscid flow. 
This may be caused by the viscous flow itself since we are 
taking into account the shear force in the boundary layer near 
the surface. 

For the afterbodies from the small perturbation solution 
Gemteour (SPS1.0, SPSi.1 and SPS1.2), in inviscid flow, the 
shocks are formed at the contour surface. The location of the 
shock depends on the specific afterbody contour and the 
approaching Mach number; at the higher Mach numbers the shock 
appears a bit further downstream as shown in Figures (16) 
through (21). As shown in Figure (22) and Figure (23), this is 
apparently a result of viscosity; it shows the boundary layer 
by the increment of Mach number away form the surface. Weak 
shocks were formed further downstream compared to inviscid 
flow. Flow separation occurs in the starting contour region 
and is followed by circulating flow in the base region. 

For the conical afterbody, the approaching free stream 
Mach number increases following the mid-section surface, then 
a flow expansion occurs at the turning angle region, until the 
maximum local Mach number is reached. A weak shock is formed 
at the end of the boattail region as shown in Figure (24), 


while for the viscous flow, the weak shock develops away from 
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boattail surface due to the boundary layer and the separation 
of the flow. The circulating flow in the base region is-more 
Significant than in the inviscid flow as depicted in Figure 
(25). In addition, from Table 2, for viscous flow the meseunum 
local Mach number for the conical afterbody is higher than for 
the contoured afterbody. 

For each afterbody the pressure drag coefiicreoneaaes 
versus free stream Mach number (M,) for inviscid and viscous 
flow are plotted in Figure (26). The negative sign GEG 
due to the fact that the calculation of pressure starts from 
mid-section through the base of afterbody and ignores the 
forebody pressure. These results show that the pressure drag 
coefficient is higher for viscous flows than for inviscid 
flows for each given afterbody. This may be caused by the 
viscosity effect and the pressure distribution difference in 
the flow field. In addition, it can be seen from the Ghare, 
that the C, for the afterbody with small perturbation Ssolmeue 
contour Mach 1.10 (SPS1.1) has the lowest C, valucs (oye ee- 
entire Mach number range. Therefore, the SPS1.1 contour shows 


to be relatively the best among these afterbodies. 


b. Complete Bodies 
Complete bodies consist of an arbitrary forebody 
(ellipsoid) joined to various afterbodies such as the small 
perturbation solution contours (SPS 1.1 and SPS 1.2) andja 


conical afterbody. The approaching free stream Mach number (M 


io) 
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Pomging from 0.50 to 1.50) Starts from outer boundary with 
distance of 1 times for the coarse grid and 5 times body’s 
length for the fine grid respectively. Significant differences 
between coarse and fine grid calculations can be seen from the 
characteristic bow shock. In the coarse grid, the bow shock 
hits the outer boundary, this causes the approaching free 
stream Mach number not to be the same as in the input 
parameter as shown in Figure (27). While in the fine grid, the 
bow shock dies out before reaches the outer boundary as 
depicted in Figure (28). So, further discussion is focused on 
the fine grid exclusively. 

The flow stagnates on the nose and then follows the body 
surface until it reaches a maximum local Mach number as 
tabulated in Table 3. The critical Mach number for these 
complete bodies 1S approximately at M, = 0.70, where the 
maximum local Mach number reaches unity at the shoulder 
region. 

Figure (29) shows a typical high subsonic free stream Mach 
number (M, = 0.85) flow over the complete body SPS 1.1. The 
flow stagnates on the nose tip, then flow is accelerated 
following the forebody surface reaching a sonic line at the 
mid-way of the forebody and forming a supersonic region at the 
shoulder. Then, the flow deccelerates at the mid-section and 
accelerates again until reaches a maximum local Mach number at 
the starting region of the afterbody and a weak shock occurs 


in this region. A supersonic approaching free stream Mach 
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number, is depicted in Figures (30) and (31) for SPS 1.1 and 
SPS_1.2 respectively. A bow shock is obviously seen in front 
of the nose. The characteristic of the bow shock is more 
inclined down stream for the higher M, as shown in Figure 
(32). A subsonic region is formed between bow shock and the 
nose, then the flow accelerates along the forebody surface up 
to a supersonic region in the mid-section. The expansion flow 
occurs in the starting contour region until it reaches a 
maximum local Mach number. Then, a shock is formed in the 
contour region. Similar as in the afterbody only, the shock 
location depends upon the contour surface and M,. The shock 
location for a given contour is more downstream for higher M, 
and at the same M,, the shock location for SPS 1.1 is more 
downstream than SPS 1.2. A typical residual history of CFD 
calculation for complete body with small perturbation solution 
contour is shown in Figure (33). Convergence is obtained at 
about 500 iterations. 

For the conical afterbody, the flow characteristic is the 
same as the other complete bodies up to the mid-section 
region. The expansion flow occurs at the turning angle, then 
the flow accelerated along the conical surface and weak shock 
is formed at the edge of base. Figures (34) and (35) show the 
Mach number contour and corresponding residual calculation for 
conical afterbody at M, = 1.10. 

The pressure drag coefficient (C,-press) versus free 


stream Mach number (M,) for fine grid complete bodies are 
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plotted in Figure (36). The drag rises sharply in the high 
subsonic Mach number (M, = 0.95) and reaches a maximum (peak) 
at M, = 1.10. Then the drag decreases with a shallow curve as 
the M, increases. The decreasing shallow curve may be caused 
by the bluntness of the nose and it agrees with Shapiro 
[Ref.18] because the fineness ratio and bluntness of the nose 
of bodies of revolution are the important factors that 
contribute the drag curve at transonic and supersonic range. 
It can be seen from the graph, the drag curve of complete body 
with conical afterbody is higher than with small perturbation 
solution contour at the entire M,. Furthermore, the peak of 
the drag curve iS approximately 15% higher. The drag curve for 
peomiel and SPS 1.2 are likely to have the same trend up to M, 
= 0.95; beyond this Mach number, the drag curve for SPS 1.2 is 
slightly greater than for SPS 1.1. Therefore, the complete 
body with small perturbation solution contour afterbody Mach 


ie ie(SPS 1.1) relatively gives the lowest drag. 


oS: 


VI. CONCLUSIONS AND RECOMMENDATIONS 


The use of numerical Simulation (CFD) appears to be the 
most cost effective method to predict the aerodynamic 
performance, especially in the transonic range. 

In this research, the grid-generating program GRAPE is 
Suitable only for the geometry with smooth contour surface 
(SPS_1.1 and SPS_1.2), while the program GRIDGEN2D is used for 
the geometry of a conical afterbody with non-smooth points. 

We have shown that fora complete body model the use of a 
fine grid (77x9x120) is more reliable than a coarse grid 
(77X13x60). This is shown by the characteristic of the bow 
shock at M, > 1 as seen in Figures 27 and 28. The pressure 
drag coefficient (Cg) versus free stream Mach number (M,) 
graphs show that the small perturbation solution contour for 
Mach 1.10 (SPS_1.1) gives relatively the lowest Cl=Gnmaaaed 
models (afterbody and complete body), a decrease of 
approximately 15% of peak in the transonic range compared to 
the conical afterbody. Therefore, the best design of an axi- 
symmetric body such as missiles, projectiles and aircraft, can 
be based on the small perturbation solution contour. 

Finally, this work may be continued with the investigation 
a complete body with pointed nose geometry and the application 
of a Euler as well as Navier-Stokes flow-solvers within 


various small angle of attack. 
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Figure 10. 2D grid, conical afterbody model. 
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Figure 11. 2D grid, complete body with conical afterbody. 
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Figure 14. 3D fine-grid, complete body with conical afterbody. 
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Figure 15. 3D fine-grid, complete body 
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APPENDIX 


COMPUTER PROGRAMS 


ve 


PROGRAM KSI.F 

KKK KK KKK KKK KKK KKK KKK KK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKKKeKEKK 
* This program is developed to calculate KSI(X), K(X) and X * 
* using mumerical integration based on trapezoidal rule to * 
* solve equation (32) and plotting the output as shown in * 
* , WAS ee ae * 
kKaekekkeke ke KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKKKK KKK KKK K KKK KKK KKKKEK 


REAL M(3), X(0:401,3), A, Al, B(401,3) 
# , KSI, P, FUNC, H, DD, Q, b, - K(40272)) ae 
INTEGER YY 

OPEN (UNIT=9, FILE=' KSI’ ,STATUS=' UNKNOWN’ ) 

PRINT *,’ENTER LOWER & UPPER BOUND, # OF INTERVALS, # OF 
#DATA SET’ 

READ *, A, Al, N, N21 

100 BOs Oei=1F2 
BRINIT*, ~ENIER MAGE Nes 


READ*, M(T) 
C M(I) = 0.9+1I*0.1 
Tee (M (1) ee 0 aameeviAN 
Pee -1. 
ELSE 
P=1. 
ENDIF 
DO 20 J=1,N1 
B(J,LI) = J*A1/N1 
H = (B(J,1I)-A)/N 
AREA = 0. 
K(J,I)=(-0.0102/ (M(I) **2) ) * ( (ABS ((B(J, I) **2) - ( (ABS (1- 
# M(I) **2))**1.7574)))**(2./3.)+ (ABS (1-M(I) **2) )**1.1716) 
Cc 
DO 30 L= 1,N 
KS] A 10m eer 
DD = (KSI**2+P* (ABS (1-M(I) **2) ) **1.7574) 
me elo ee rte) wena 
FUNC = 1./(DD**(1./3.)) 
QO= 1 
ELSE 
FUNG = 1../ (( 4d DD) =~ (ieee a 
QO = -1. 
ENDIF 
AREA = AREA + H*FUNC*Q 
(Ciera ee 
30 CONTINUE 
20 CONTINUE 
16 CONTINUE 
(C 


DO 40 J= 1,N1 
PRINT 50, (Bld, 2), X\0; 1) ee ee 


78 





a8, FORMAT UX 3 (bG.3 (201075, 1X) ) 
40 CONTINUE 
C 
200 PRINT*,’TRY AGAIN?, ENTER 1 FOR YES,2 FOR DATA FILE, 
# OTHERS FOR NO’ 
READ*, YY 
They nO. 2) TEN 
GO TO 100 
ELSE 
eo. BO) 2) SHEN 
BO Le t=" 1 Ni 
While (9750) )( Bi, 1),2,1),ki7,1) ,1=1,3)} 
L1G CONTINUE 
GO TO 200 
ENDIF 
ENDIF 
END 


PROGRAM ZETA.F 
KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKEKKKKKKKKKKKKKKKKKKKK 
* This program is written to calculate ZETA(r) and r using* 


eecticat+on (29) and plotting Output an Fig.3 . * 
KKKKKKKKKKKKKKKKKKEKEKEKEKEKKKKKKEKKEKEKKE KK KK KKK KKKKKKKKKKKKKKKKKKEK 


REAL M(3), ZETA(0:14,3), R(14) 
OPEN (UNIT=9, FILE=’ ZZ’ , STATUS=’ UNKNOWN’ ) 


Cc 
DO 10 I= 1,3 

cS M(I) = 0.7+1*0.1 
PRINT*, ’ENTER MACH NO.’ 
READ*, M(TL) 
DO WwZ0 J= 1,14 
Zee. .) = 2000.0 
PRINT*, ‘ENTER R VALUES’ 
READ*, R(d) 

a! Rig) = 0. Pad 


ZETA(J,1I) = (4/(R(dJd) **2)) + (ABS (1- (M(I) **2) ) 
# *(R(J) **2.8284) ) 
# +((1-M(1I) **2)*(R(J) **7.657) /50.63 
IF (ZETA(J,1) .GE.ZETA(J-1,1I) ZETA(J,1)=ZETA(J-1,1I1) 
20 CONTINUE 
me CONTINUE 


PRINT 30, M 
30 FORMAT (11X,5F10.4) 

DO 40 J= 1,14 

WRITE (9750), (RR (i )5 (ZETA), 1), 2=1,4) 
50 FORMAT (4X, 3 (2F10.4, 2X) 
40 CONTINUE 

END 


fas, 


PROGRAM D2D3.F 
KKKKKKKK KKK KKKEKEKKEKKKKKEKKKE KKK KKK KKK KK KKK KK KKK KKK KR KKK KKK K 
* This program will read data files output Gridgen2d * 


* from 2D to 3D for further processing with rotategr.f * 
KEK KKK KKEKEKKKKEKEKKKKRKEKKKKKEKKKRKE KK KKK KRKRKRKKKKRKKKRKKKKKRKKKRKRKRKRKKK KEK 


REAM X(77,60N Y( 79) GON ae oD 

READ (12,*) IDIM,JDIM 

READ (12,*) ((X(I,J),I=1,IDIM) ,J=1,JDIM), 
# ((Y(L,d) , l=1, 1DIM) ,J=—, JDm) 
Zi, J oO 

WRITE(14,*) IDIM,JDIM,1 

WRITE(14,*) ((X(1I,Jd),I=1,IDIM) ,J=1,dJDIM), 
# ((0.0,I=1,IDIM) ,J=1,JDIM), 
# (CY (LJ), 1=1) 1D EM ae — dee 
STOP 

END 


PROGRAM HALF.F 
KK KKK KKK KKK KK KKK KK KEK KEKE KKK KKK KKK KKK KEK KEKE KEKE KEKE KKEKKEKKKEKKEKE 
* This program will write out half of a grid file output * 


* from Grape for further processing with rotategr.f * 
KKK KKK KKK KKK KKK KKK KKK KEKE KEKE KKK KEKE KKKEKEKEKKKEKKKKKKKEKKEKKEKKKEKE 


DIMENSION X(2007 200) Y (2007200) 

READ (30,*) IDIM,JDIM 

READ (30),*)  ((( 1,0) 3-1 oo ae 
J=1 





# ((Y(I,J),I=1,IDIM) , J=1,JDIM) 
ts 
IDMD2=IDIM/2 + 1 | 
Cc i} 
DO 10 J=iegbiM " 
Y{(IDEMe) go) =80.0 iF 
10 CONTINUE | 
Es 
WRITE (31,*) IDMD2,JDIM,1 | 
WRITE (31,*) ((X(I,J),I=IDMD2) ,J=1,JDIM), 
# ((0.0, I=IDMD2) ,J=1,JDIM), 
# ((Y (I,J), I=IDMD2) , J=1,JDIM) | 
STOP 
END 


80 


Sl 


PROGRAM ROTATEGR.F 


kKaekekeKkkeke Keke KKK KK KKK KKK KK KKK KKK KKK KKK KKK KK KKK KK KKK KKK KKK KKK KK KK 


* This program will create 3D grid (body of revolution) 
* by rotating 3D plane grid for further processing with 


* Flow Solver 
kek kk kkkeKkeKeKe KKK KK KKK KKK KKK KKK KK KKK Ke KKK K KKK KKK KKK Ke KKK KKK KKK KEKEKK 


oe 


ial 


DIMBNSTON X(200;,200)5.. Y(200;.200) 7 2(200;,200) 
BiIMeNslon xx (200,100,200)4 YY(200, 100,200), 
Hoa i200, 100.200) 


CHARACTER*30 FNI 
CHARACTER*30 FNO 


PRINT*, ‘INPUT GRID’ 
READ (*,21) FNI 


PRINT*, ‘SG = 0, MG = 1’, ‘’',UNF=1, FORM=2' 
READ (*,*) IGRI 

PRINT*, ‘IFORI = 1 UNF, IFORI = 2 FORM’ 
READ (*,*) IFORI 

FORMAT (A) 

Teen PORT-EO.1) THEN 

REWIND 1 


OPEN (1, FILE=FNI, FORM=’ UNFORMATTED’ ) 
IF (IGRI.EQ.1) READ (1) MGR 
READ (1) “Tl, Ji, K1 
Ren sole (CLL wd) abate) Hla), 
CY (rd) a al 
CUZ eG tal) p= at) 
CLOSE (1) 
ENDIF 
IF (IFORI.EQ.2) THEN 
REWIND 2 
OPEN (2, FILE=FNI, FORM=’ FORMATTED’ ) 
IF (IGRI.EQ.1) READ(2,*) MGR 
READ(2 ~*) “1L,J51,41 
READ (27%) Cixi ad) ha P= 1), 
Cy) ) rope ape ea pel lee ) ’ 
CZ oe ee a ieee a — 
GHOSE \(Z) 
ENDIF 


PI = 4.*ATAN(1.) 

PRINT, NO OF PLANE -IN J DIRECTION 7?’ 
READ(5,*) JM 

DTH = (180./(JM-1)) * (PI/180) 

iO @ Mie ie as ees a 

DOL li d= 0 


K=J 

ROG, 24.0) (=) Kt) 
Wt 2c = 0.0 
Zia, op), = ZL, K) 
CONTINUE 


ail 


* 
* 
* 


Zo 


20 


# 
# 


TM=I1 
KM=J1 

DO 20 J=3,JM+1 
DO 20 I=1,IM 

DO 20 K=1, KM 

RX (Tdi KR) =x (eK) 


Abie. (te) 3h) DEL 

YY(i jd; hk) = Sathya) 
ZZ4i, 307) = COs th) += 2a 
CONTINUE 


DGe20 2-1, 1M 
DO 30 K=1,KM 


a 

XX (Eph) = eee +2) KX) 
VYY¥(r Jaki =- YY (Evers, 
Zit, J pie = 4A lee, Kh) 
J=JM+2 

PORE ig) =e (eG 2 
YY Ul wah) =e (ee) 
LAE , Kh) See be) 
CONTINUE 

JM=JM+2 


PRINT*, ‘Output filename =’ 
READ (5,21) FNO 


REWIND 3 

OPEN( 3, FILE=FNO, FORM=’UNFORMATTED'’ ) 

WRITE (3) IM,JM, KM 

WRITE(3) (((X(1I,d,K), I=1 40M) g=s990M) | Kee 
(( OYOCE , dake, Tai, 1M) 02 oie) 
(0 Za lene) ay t= 2 ee = ee ee 

Chess (3) 

STOP 

END 


82 























PROGRAM READX.F 


KKEKKKKEKKKKEKEKEKEKKKKKKKEKKEKKKKKKEKKEKKKKK KKK KK KKK KKK K KK KKK KKEKKEKK 


* This program will read 3D grid file and convert from * 
* formatted into unformatted file (grid.for to grid.in) * 
* for further processing with Overflow. * 


kKkekekkKeKkKeKeKKKeKKKKKKKEK KKK KKK KKK KKKKKKKEKKKKKKKKKEK KKK KKK KKK K KEK EK 


DIMBNSEON % (77,13, 60)ry (97. 13.460), 2i77,43660) 
OPEN (UNIT=12, FILE=’GRID.FOR’ , STATUS=' UNKNOWN’ ) 
OPEN (UNIT=14, FILE=’GRID.IN’ ,STATUS=’ NEW’, 
#EFORM=’ UNFORMATTED’ ) 

READ (12,*) IDIM,JDIM, KDIM 

READ (12,20) (((X(I,J,K), I=1,IDIM) 
# (((¥(I,J,K), I=1, IDIM) 


, , JDIM) ,K=1,KDIM) , 
( ,J= 
# (((Z(I,d5,K), I=1,IDIM) ,J 


=1 
bor ara ee 
=1,JDIM) , K=1, KDIM) 


WRITE (14) IDIM,JDIM, KDIM 

WRITEUPS) AO, Ke 2=1 7 Dim). J=1, IDEM), K=1,.KDIM) , 

# COG (eK) lh, IDEM) (oOo =), JD IM), Ki KD IM). 

# CA td, we 2sl, LeiM) j,o=-1,JDIM), K=1, KDIM) 
ZU FORMAT (6F15.6) 

SLOP 

END 


PROGRAM READQ.F 


kKekekkekeKeKkeKkeKeKeKKeKeKeKeKK KKK KKKKKKKKKKEKKKKKKEKKKEKKEKKEKKKEKKKKKKKKKKKKK 


* This program will read q.save file output from * 
* Overflow and convert from unformatted into formatted * 
* file for further processing with Plot3d. * 


kekkkekKeKkKkeKeKeKKeKeKKKeKKKKKKKK KKK KKK KKK KKK KKK KKKKKKKKKKKKKKKEKKEK 


DEMENSION ©0(77,13,60, 5) 

OPEN (UNI T=. FIRE= O.SAVE’ ,STATUS=’ OLD’, 

#FORM=’ UNFORMATTED’ ) 

OPEN (20, FILE=’Q.FORM’ ,STATUS=' NEW’ , FORM=’ FORMATTED’ ) 
READ (1) NI,NJ,NK 

READ (1) FSMACH, ALPHA, RE, TIME 

READ (1) A VOT NA) peer NL) o=1, Ne) UK) NK N= 1S) 


WRITE(20,*) NI,NJ,NK 
WRITE(20,*) FSMACH, ALPHA,RE, TIME 
WRimds (20 2* CC (Ol Lg KNX) 


fol No) ke NK) | XH 1; 5) 
PRINT*, NI,NJ,NK 

SOL 

END 
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PROGRAM OVERFLOW. IN 
Kae K KKK KK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK K 
* Input file for Overflow to run either EULER or NAVIER- * 
* STOKES Flow-seiver. « 
Keke KK KK KK KKK KKK KKK KKKKKKKKKKKKK KKK KKK KKK KKK KKK KKK KKK KR KKK KKKK 
SGLOBAL 
CHIMRA=.F., NSTEP=100, RESTART=.F., NSAVE=100, NOT=0, 
SEND 
SFLOINP 
ALPHA=180.0, FSMACH=0.95, REY= 6.00E6, TINF-= 5202 
SEND 
SGRDNAM 
NAME= ‘'AXI-SYMMETRIC BODY WITH AFT-CONE, 77x13x60 
GREED, 
SEND 
SMETPRM 
TRHS = 0, ILBS = 92) 22) lSS a, 
SEND 
STIMACU 
DT = 0.1, ITIME = 1, TFOSO = 1.00, CELMIN — 2s 
SEND 
SSMOACU 
ISPECJ 
ISPECK 
ISPECL 
SMOO = 1. 
SEND 
SV lo LN 
VISCJI 
NTURB 
IESE Ae 
POT 
JTLS 
JTLE 
KTLS 
KTLE 
LTLS = 
LTLE = 
TLPAR1= 
SEND 
SHS OIIN NS 
NBC = 
ie iy P 
16) =) BAR! 
JEBCS. >= oe de i 
JBCE 
KBCS 
KBCE 
LBCS = Ls 
LBCE = S10) fume Sy (S10) 
SEND 


we, 
HH 
Com ea) 
NN 
AG 
| 
ONNN 


00) ~DES au 
.00, DIS4K 
.00, DIS4L 
Sia ys 


meear 
oO 
OZ, 


I 
Lai (i an Le 


Nol 

ONNN 
O 
IH 
ep) 
N 
i 
I 


. VISCK =p we Vacer 


nou uw 
lI 
rr 
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I Il 
e t- ~] e 
WOOF WR QR We O A] 
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Nae 
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